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Abstract. We study a nonlocal version of the one-phase Stefan problem 
which develops mushy regions, even if they were not present initially, a model 
which can be of interest at the mesoscopic scale. The equation involves a 
convolution with a compactly supported kernel. The created mushy regions 
have the size of the support of this kernel. If the kernel is suitably rescaled, 
such regions disappear and the solution converges to the solution of the usual 
local version of the one-phase Stefan problem. We prove that the model is 
well posed, and give several qualitative properties. In particular, the long- 
time behavior is identified by means of a nonlocal mesa solving an obstacle 
problem. 



The aim of this paper is to study the following nonlocal version of the one-phase 
Stefan problem posed in x (0, oo), 



which presents interesting features from the physical point of view. The function 
J is assumed to be continuous, compactly supported, radially symmetric and with 
/uN J = 1- We denote by Rj the radius of the support of J. For some results we 
will also assume that J is nonincreasing in the radial variable. 

The local model - The well-known usual local Stefan problem is a mathematical 
model that describes the phenomenon of phase transition, for example between 
water and ice, [H], (35]. Its history goes back to Lame and Clapeyron and, 
afterwards, Stefan [23'. The one-phase Stefan problem corresponds to the simplified 
case in which the temperature of the ice phase is supposed to be maintained at the 
value where the phase transition occurs, say 0°C. The thermodynamical state of 
the system is characterized by two state variables, temperature v and enthalpy u. 
Conservation of energy implies that they satisfy, in the absence of heat sources or 
sinks, the evolution equation 



where the density p > and the thermal conductivity n > are assumed to be 
constant. 
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On the other hand, there is a constitutive equation relating v to u, given in an 
ideal situation by the formula 

(1.2) v = c-\u-L)^, 

where c > 0, which we assume to be constant, is called the specific heat (the amount 
of energy needed to increase in one unit the temperature of a mass unit of water) 
and L > is the latent heat (the amount of energy needed to transform a mass unit 
of ice into water). All the parameters can be set to one with a change of units, and 
we arrive to 

(1.3) dtu^l^v, v = {u-l)+. 

In contrast with the standard heat equation, this problem has finite speed of prop- 
agation: if the initial data are compactly supported, the same is true for u (and 
hence for v) for any later time. This is one of the main features of the model, and 
gives rise to the existence of two free boundaries, one for u, the set d{u > 0}, and 
one for v, d{v > 0}. 

On the ice we have u = 0, while u ^ 1 on the liquid phase. The points where 
< u < 1 correspond to the mushy region, where we have neither ice nor water, 
but something in an intermediate state between solid and liquid. The temperature 
in this zone is 0, but not its enthalpy. 

There is a major drawback in the model: either u{x, i) ^ 1 or u{x, t) = f{x) [15j . 
This means, on one hand, that a point can belong to the mushy region only if it 
belonged to it at the initial time. On the other hand, there is no evolution of the 
enthalpy inside the mushy region. Both things are unsatisfactory from the physical 
point of view. Indeed, once the structure of ice starts to break, it should take some 
time until it melts completely. The region where this occurs may be small, but it 
should be noticed at some intermediate (mesoscopic) scale between the microscopic 
and macroscopic ones. 

The nonlocal model - Bearing the above discussion in mind, we consider the 
nonlocal diffusion version (jl.ip of the Stefan problem. As we shall see, this model 
also has finite speed of propagation, and hence free boundaries. Moreover, in sharp 
contrast with the local problem, it develops mushy regions, even if they were not 
present initially. The size of these mushy regions is given by the size of the support 
of the convolution kernel J . Finally, if the initial data are continuous, the solution 
remains continuous for all times. 

It is possible to rescale the kernel in such a way that solutions of the nonlocal 
model converge to solutions of the local version ()1.3p . In this scaling, the support 
of the kernel shrinks to a point, and hence the mushy regions disappear. Thus, the 
local model can be viewed as a limit problem when going from the intermediate 
scale to the macroscopic one. 

What is the point in considering this particular diffusion operator, Cv :— J^v—v? 
With this choice, the evolution of u at a certain point is governed by a balance 
between the value of the temperature at this point and a certain average of this 
physical magnitude in a fixed neighborhood. This is a way to take into account 
possible middle-range interactions between water and ice. In the local model, the 
region where the average is taken shrinks to a point, and the evolution of the 
enthalpy is governed by the Laplacian of the temperature. 
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Choice of the kernel - One could use other non-local operators, as fractional- 
Laplacian type operators of the form Kv —{~/^Yv. However, we will stick 
to operators C involving kernels which are both compactly supported and non- 
singular, though most (if not all) of our results should be valid for singular (yet 
compactly supported) kernels. Why are we being so restrictive? There are two 
main reasons: 

{i) Using a non-compactly supported kernel means that "infinite-range" inter- 
actions are not considered to be negligible. This would imply for instance 
that the mushy regions which are automatically created are instantaneously 
spread throughout the space (see Section 1531 with Rj = oo), something 
which is not satisfactory. 

{ii) Kernels which have a singularity at the origin have the big advantage of 
implying regularization properties ([6] and |11|). It is challenging from a 
mathematical point of view to try to handle situations where no regularizing 
effect can occur. 

Another way to think about this absence of regularization is to under- 
stand it as a lack of compactness of the inverse of £, which is the source of 
several problems. Our approach, which consists in controlling the supports 
of solutions, allows us to bypass such problems, up to a certain point. 

About the initial data - We would like to consider typical situations in which 
initially there is ice with zero enthalpy everywhere except at some places where 
there is water at some positive temperature. In this case, we shall see that a 
mushy region will develop as ice melts from an initial configuration in which no 
mushy region was present. Thus, we have to consider initial data which are in 
principle not continuous, but only integrable. This will lead to solutions that are 
not continuous. Therefore, we have to handle the various qualitative properties 
concerning the supports and mushy regions not in the usual (continuous) sense, 
but in the sense of distributions, which requires some effort. 

If the initial data are continuous and bounded, we can construct a solution with 
these two properties. Some of our techniques can be greatly simplified for such 
solutions, and can be later applied to more general data by approximation. Hence 
we will devote some time to discuss the basic theory for solutions in this class. 

As expected, when the initial data are at the same time bounded, continuous 
and integrable, the two concepts of solution mentioned above coincide. 

In the one-phase Stefan problem, the temperature, and hence the enthalpy, are 
assumed to be nonnegative. Therefore, in principle we should only consider non- 
negative initial data /. Indeed, if the initial enthalpy had sign changes, we would 
be dealing with a two-phase Stefan problem, and the relation between the temper- 
ature and the enthalpy would not be given by (|1.2p . In particular, the temperature 
should be negative for negative values of the enthalpy. Therefore, the relevant 
model would be different. 

Though the problem might not have any physical meaning for functions / that 
have sign changes, some of our results are still true for such general initial data. 
Whenever we know that this is the case, we will state the corresponding theorem 
for a general /, and will add the restriction of nonnegativity only if required. 
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Long-time behavior - Another important aspect of evolution equations is the 
long time behavior of solutions. For a certain class of initial data (which does not 
include all non-negative functions in L^(K^) if = 1, 2), we prove that the solution 
converges as time goes to infinity to a solution of a nonlocal mesa problem which 
is not the "classical" (local) mesa problem. The identification of the limit (the 
mesa) can be done by solving an obstacle problem. More precisely, the asymptotic 
behavior of the solution u is given by the projection 

Vf :^ f + J * w — w 

where w solves an obstacle problem detailed in Section [5] 

The projection operator V is L^-contractive in the set S of admissible functions 
for which the above convergence result is valid, CoroUarv 15.51 On the other hand, 
S is dense in the set of L^(M^) functions which are non-negative a.e. Therefore, 
V can be extended by continuity to the latter class of functions. This will allow us 
to identify the large time limit of solutions of for the wider class of integrable 
and nonnegative initial data. 

Abstract setting - Equations like (II. ip and (|1.3p can be embedded in the ab- 
stract setting of semi-group theory for equations of the form 

where 7 is a maximal monotone graph, and ^ is a linear (bounded or unbounded) 
TO-accretive operator, see [T5^ and also [1], [5]. This theory provides existence and 
uniqueness for our model, though not the main qualitative properties considered 
here. As for existence and uniqueness, in the special case that we are studying, 7 
is the inverse of the Lipschitz graph F : s n> (s — 1) + , and many arguments can be 
written at a lower cost. 

Notation - Given a set we define: 

• BC{n) ^{ipe C{n) : ip bounded in R^}; 

• Cc(ri) = {if E C(ri) : (f compactly supported}; 

• C^(r2) ^ {^p e C°°{Q) : If compactly supported}. 
We also denote 

• L^(M^) ^{ipe Li(M^) -.ip^O a.e.}; 

• BC+(M^) BC(M^) -.ifi^O}; 

• Co(M^) = {fe C(M^) : ^ as |a;| -> 00}. 

The elements belonging to L^((0, T); L^(M^)) will sometimes be viewed as 
elements of L-^(R'^ x (0,T)). In such cases we denote tp{x,t) = ip{t){x). 

Organization of the paper - In Section 2 we derive the general theory of 
the model both for integrable initial data and for continuous and bounded initial 
data. Section 3 is devoted to the study of mushy regions and free boundaries. 
Convergence to the local Stefan problem and disappearance of the mushy regions 
in the macroscopic scale are done in Section 4. In Section 5 we study the large time 
behavior of solutions. Numerical experiments and illustrations of the qualitative 
properties of the model are collected in Section 6. We devote the last section of the 
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paper to establish some conclusions and make some comments on the model. In 
order to prove some of our results, we have needed to improve slightly the existing 
results on the asymptotic behavior of solutions to the non-local heat equation. Such 
improvements are proved in an appendix. 

2. Basic theory of the model 

We will develop here the basic theory for the two concepts of solution mentioned 
in the introduction. 

2.1. theory. We start with the theory for integrable initial data. In this case 
the solution is regarded as a continuous curve in L^(M^). 

Definition 2.1. Let / e Li(R^). An L^-solution of (HU) is a function u £ 
C([0, oo); L^(]R^)) such that ()1.1|) holds in the sense of distributions, or equiva- 
lently, if for every t > 0, u(t) e L^{R^) and 

(2.1) u{t) ^ f + [ {J *v{s) -v{s))ds, v^{u-l)+ a.e. 

Remark. If u is an L^-solution, then u € L'^{R^ x [0, T]) for all T > 0. Hence, ([rT|) 
holds, not only in the sense of distributions, but also a.e., and u is said to be a 
strong solution. Moreover, since v — {u — 1)+ G C([0, oo); L-'^(K-'^)), we also have 
u G C^([0, oo); L^(R^)), and the equation holds a.e. in x for all t ^ 0. 

Theorem 2.2. For any f e L^(R^), there exists a unique 'L} -solution of ()1.1|) . 

Proof. Let Bt^ be the Banach space consisting of the functions u e C([0, t^] ; L^ (M^)) 
endowed with the norm 

|||m||| = inax ||m(<)||li(r«). 
We define the operator T : Bt^ Bt^ through 

{Tfu){t) ^f+ f\j*{u~ l) + {s) ~{u- 1)+(,s)) ds. 
Jo 

This operator turns out to be contractive if to is small enough. Indeed, 

< / r (\J * ((^ - 1)+ - ('A - 1)+)(s)| + - 1)+ - (V^- - 1)+|(s)) ds 
Jk« Jo ^ ' 

< ^ (ll J||li(r«) + 1)) II((<P - 1)+ - (V^ - 1)+)(s)||li(r«) ds. 

Hence 

IIIT/^ - T/^lll 2io max ||((^ - 1)+ - (^ - l)+)(0 1|lmk«) 

Thus, T is a contraction if tg < 1/2. Existence and uniqueness in the time interval 
[0, to] now follow easily, using Banach's fixed point theorem. Since the length of 
the existence and uniqueness time interval does not depend on the initial data, we 
may iterate the argument to extend the result to all positive times. □ 
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The energy of the L^-solutions is constant in time. 
Theorem 2.3. Let f e L^{R^). The -solution u to dTH]) satisfies 

"(^) — I f fo''' every t > 0. 



Proof. Since u{t) e L^(R^) for any t, integration of the equation (j2.ip in space 
yields, thanks to Fubini's Theorem: 

uit)^ f f+f(f J*v{s)-( vis))ds^ f f. 

□ 

Our next aim is to derive an -contraction property for L^-solutions. In order 
to obtain it, we need first to approximate the graph r(s) = (s — 1)+ by a sequence 
of strictly monotone graphs r„(s) such that: 

(i) there is a constant L independent of n such that |r„(s) — r„(i)| ^ L\s — 1\, 
for aU n e N and s,t^O; 

(ii) for all n G N, r„(0) = and r„ is strictly increasing on [0, oo); 

(iii) r„(s) ^ s for ah n £ N and s ^ 0; 

(iv) r„ — ?• r as n ^ oo uniformly in [0, oo); 

(take for instance r„(s) = s/(n + 1) for ^ s ^ (n + l)/n, and r„(.s) = s — 1 for 
s > (n + l)/n). 

Since r„ is Lipschitz, for any / e L^(R^) and any n G N there exists a unique 
L^-solution Un G C([0, oo); L^(R^)) of the approximate problem 

(2.2) dtUn = J * r„(u„) - r„(M„) 

with initial data w„(0) = /. The proof is just like that of Theorem 12.21 Moreover, 
r„(it„) e C([0,oo);Li(R^)), and, hence, w„ G CH[0, oo); Li(R^)). Conservation 
of energy also holds. 

The -contraction for our original problem will follow from an analogous result 
for the approximate problems. 

Lemma 2.4. Let Un,i and Un,2 be two -solutions of ^2.2\ with initial data /i, /2 £ 
Li(M^). Then, 

(2.3) / {un.i - Una) Ai) I (/1-/2), for every f^i Q. 

Proof. We subtract the equations for Un,i and u„^2 and multiply by i>m„2}- 
Since u„,i - it„,2 G Ci([0, 00); Li(]R^)), then 

dt{Un,l - Uri,2)l{«„,i>u„,2} = dt{Un,l - Wn,2)-|-- 

On the other hand, since ^ l{ti„ i>u„ 2} ^ 1' have 

J * (r„(M„,i) - r„(u„,2))i{u„,i>u„,2} ^ J * (r„(u„4) - r„(u„^2))+- 

Finally, since r„ is strictly monotone, l{„„,i>„„_2} = l{r„(«„,i)>r„(«„,2)}- Thus, 

(r„(w„,i) - r„(u„,2))i{«„_i>«„,2} = (r„(M„,i) - r„(M„_2))+- 
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We end up with 

dt{uns - Mn,2)+ ^ J* (r„(u„,i) - r„(u„_2))+ " (r„(u„,i) - r„(w„_2))+- 

Integrating in space, and using Fubini's Theorem, which can be appUed, since 
(r„K,i(i)) - r„(w„,2(t)))+ G Li(R^), we get 

dt / [Un.l - Un,2)+(<) ^ 0. 

□ 

Remark. It follows immediately from (j2.3p that 

||(u„,i - W„,2)(i)||Li(R") ^ ll/l ~ /2|1l1(R")- 

Corollary 2.5. Let ui and U2 he two -solutions of (jl.ip wit/i initial data /i, /2 G 
Li(R^). Then, for every t > 0, 



/( 



(2.4) / {u,-U2)+{t)^ / (/i-/2)+. 



Proof. The idea is to pass to the limit in n in the contraction property (j2.3p for 
the approximate problems (j2.2p . Hence, the first step is to prove that any solution 
u of (ll.ip is the limit of solutions it„ to (|2.2p . 

Let w be an open set whose closure is contained inM^x(0,cx3), w CC R"x(0,oo). 
By the conservation of energy, ||Mn(OllLi(R«) = II/IIli(r")- Hence {w„} is uniformly 
bounded in L^(a;). Therefore, in order to apply Frechet-Kolmogorov's compactness 
criterium, it is enough to control 

/ = jj \un{x + h,t + s) — Un{x,t)\dx dt 

for h and s small enough (how small not depending on n). 
On one hand, thanks to the L^-contraction property, 

/ \un{x + h,t + s) — Un{x,t + s)\dxdt 

^R" 

^11 \f{x + h)- f{x)\dxdt<:Ton{l) 

Jo JR" 



(2.5) 



as /i — > uniformly in s and n. On the other hand, using the regularity in time, 
then Fubini's Theorem, and finally the sublinearity of r„ and the L^-contraction 
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property, we get 

\un{x,t + s) — M„(x,<)| Axdt 



-'R™ 



(2.6) 



f f r 

^ / \dtUn\{x,T)dTdxdt 

JO JR™ Jt 

I J * r„(u„) - r„(u„)|(x,T) dxdr dt 



t+s 



= (R") + l)llr„(u„(T))||Li(B")dTdt 



^ (IkIIl°=(r«) + i)II/IIli(r")«^- 

Taking T such that w C x (0,r), and using the estimates (|2?5l) and (|2^ we 
get the required control. 

Summarizing, along a subsequence (still noted m„ — )■ u in Lj'Q^(M^ x (0, oo)) 
for some function u. Moreover: (i) since the sequence {unit)} is uniformly bounded 
in L^(K^), we deduce from Fatou's lemma that for almost every t > 0, u{t) e 
L^(R''^); (ii) using that the nonlinearities r„ are uniformly Lipschitz, and their 
uniform convergence, we get that r„(-u„) ^{u) in Ljq^(R^ x (0, oo)); (iii) as a 
consequence, since J is compactly supported, J * r„(u„) -> J * T{u) in Lj'Q^(R^ x 
(0, oo)). All this is enough to pass to the limit in the integrated version of (|2.2p . 

Un{t) = f + {J * r„(u„(s)) - r„(u„(s))) ds, 

for almost every t > 0. If we extend M(t) to all t > by continuity, so that it 
belongs to the space C^([0, oo); L^(M^)), we get that u is the L^-solution to (II. ip 
with initial data /, i.e., u = u. As a consequence, convergence is not restricted to 
a subsequence. 

Now we turn to the contraction property. Let ui, U2 be the L^-solutions with 
initial data /i and /2 respectively. We approximate them by the above procedure, 
which yields sequences {un,i}, i = 1,2, such that Un.i Ui in Lj'Q^(M^ x (0, oo)) 
(and hence a.e.). The approximations satisfy (12. 3p . Using Fatou's lemma to pass 
to the limit in this last inequality we get that (12. 4p holds for almost every t ^ 0. 
Finally, since the solutions are in C([0, oo); L^(R^)), we deduce that this inequality 
holds for any t ^ 0. □ 



Remark. Equation (j2.4l) implies a comparison principle. In particular, if / 5^ 0, we 
conclude that u ^ 0. Hence, u is truly a one-phase solution. 

The temperature turns out to be subcaloric. 

Lemma 2.6. Let f G L^(R^) and u the corresponding -solution. Then the 
temperature u = (w — 1)+ satisfies Vt ^ J * v — v in the sense of distributions and 
a.e. m x (0, oo). 



Proof. Since u G Ci([0,oo);Li(R^)), 

Vt = utt{u>i} a.e. 
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In the set {u ^ 1} we have ut — J * {u — 1)+ ^ and vt ~ 0, whereas in the set 
{u > 1} we have vt — ut a.e.. In both cases we obtain vt^ut — J*v — v a.e., and 
in the sense of distributions since these are locaUy integrable functions. □ 

This property ahows to estimate the size of the solution in terms of the L°°-norm 
of the initial data. 

Lemma 2.7. Let f e L\M.'^) n L°°(R^). Then the L^-solution u of (HH) sat- 
isfies ||w(t)||Loo(Rjv) ^ ||/||l~(r") for any t > 0. Moreover, limsupj^^ -(/(t) ^ 1 
a.e. in M^. 

Proof. The result is obvious if ||/||l°°(r") ^ Ij since in this case u{t) — f for any 
t > 0. So let us assume that ||/||l°°(r") > 1- Since v is subcaloric and locally 
integrable we may use [8, Proposition 3.1] (with a compactly supported kernel), to 
obtain that 

^ ||w(i)||L~(R") ^ IK/ - 1)+IIl~(r") = II/IIl~(r«) - 1- 
Therefore, ||Loo(RiV) ^1 + ||w(i)llL~(R") =^ II/IIl~(r")- 

We can also compare v with the solution w of the following problem: 

wt^J*w^w, w{0) ^{f- 1)+ e L\R^) n L°°(M^). 

Using Theorem lA.ll in the Appendix, we obtain that the solution v goes to zero 
asymptotically like rf^^/^, so that {u — 1)+ — )■ almost everywhere, which implies 
the result. □ 

2.2. BC theory. We now develop a theory in the class of continuous and bounded 
functions whenever the initial data / belong to that class. 

Definition 2.8. Let / £ BC(M^). The function m is a BC-solution of (dU if 
u e BC(M^ X [0, T]) for ah T G (0, oo) and 

u{x,t) ~ f{x) + / [j * v{x, s) — v{x, s)) ds, w = (u — 1)+ 
Jo 

for all x e and t e [0, oo). 

Notice that if u is a BC-solution, then ut is continuous. Hence equation (|l.ll) is 
satisfied for all x and t, and w is a classical solution. 

Theorem 2.9. For any f e BC(M^) there exists a unique BC-solution of (|l.ip . 

Proof. As in the case of L^-solutions, existence and uniqueness follow from a fixed- 
point argument. 

We start by proving existence and uniqueness in a small time interval [0, to]- We 
define the operator T : BC(R^ x [0,to]) -> BC(R^ x [0,to]) through 

iTfu){x, t) = f{x) + [ ( J * (li - l)+(a;, s)-{u- l)+{x, s)) ds. 
Jo 
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This operator is contractive if < 1/2, which implies the local existence and 
uniqueness result. Indeed, a similar computation to that of the proof of Theorem l2.2l 
yields 

\Tf(p — Tfiplixjt) 2to max — M, t E \0,to]- 

R"x[0,to] 

By iteration, taking as initial data u{x,to) S BC(M^), we obtain existence and 
uniqueness for [0, 2io] and hence for all times. □ 

The BC-solutions depend continuously on the initial data. 

Lemma 2.10. Let ui and U2 he the BC-solutions with initial data fi, /2 G BC(M^) 
respectively . Then, for all T G (0, oo) there exists a constant C = C{T) such that 

max |ui - U2\{x,t) C(r)max|/i - /a], t e [0,7]. 



Proof. Since Ui is a fixed point of the operator 7/;, we have (see the proof of 
Theorem [ 



\ui - U2\{x,t) f^\fi - f2\{x) +2tQ max |ui-U2|, x G M^, t e [0, tg]- 

R"x[0,to] 

Taking t^ = 1/4, we get 

max lui — M2 1 ^ 2 max I f i — fa I , 

R"x [0,1/4] K" 

from where the result follows by iteration, with a constant C(T) = 2^-^. □ 

We also have a control of the size of the solutions in terms of the initial data. 
The proof is identical to the one for L^-solutions. 

Lemma 2.11. Let u be the BC-solution u of (jl.ip with initial data f E BC(K^), 
and let v be the corresponding temperature. Then, v is subcaloric, and \\u{-, t) ||l°°(r") ^ 
II/IIl~(R") for any t > 0. 



3. Free boundaries and mushy regions 



In the sequel, unless we say explicitly something different, we will be dealing with 
L^-solutions. Since the functions we are handling are in general not continuous in 
the space variable, their positivity sets have to be considered in the distributional 
sense. To be precise, for any locally integrable and nonnegative function g in 
M^, we can consider the distribution Tg associated to the function g. Then the 
distributional support of g, suppj5/(5) is defined as the support of Tg-. 

suppp, {g) '■— \ O, where O C is the biggest open set such that Tg\o = 0. 

In the case of nonnegative functions g, this means that x G suppjj, (5) if and only if 

G C,°°(R^), ^ ^ and ^{x) > ^ / 9{y)v{y)dy > 0. 

JR« 

If g is continuous, then the support of g is nothing but the usual closure of the 
positivity set, suppx,i{g) = {g > 0}. 
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3.1. Existence of free boundaries. We first prove that tiie solution does not 
move far away from tlie support of v. 

Lemma 3.1. Let f Then, 

suppp,(ut(t)) C suppp,(w(i)) + Br, for any t ^ 0. 

Proof. Recall first that the equation holds down to i = so that we may consider 
here t ^ (and not only t > 0). Let e C^(A'=), where A = suppp,(w(i)) + Brj. 
Notice that the support of J * u (which is a continuous function) lies inside A, so 
that 



(J * v)ip = 0. 

r 

Similarly, the supports of v and ip do not intersect, so that 

UtV — (J * v)ip — v(p — 0, 
which means that the support of ut is contained in A. □ 

As a direct consequence, we get the finite speed of propagation property. 

Theorem 3.2. Let f G Ll'^(M^) and compactly supported. Then, for any t > 0, 
the solution u{t) and the corresponding temperature v{t) are compactly supported. 

Proof. Estimate of the support of v. Notice first that 
{J*{u- l)+){x,t) < II J||l~(k")II(m - 1)+||li(r«) ^ ||>/||l~(r")II/IIli(k") := Co, 

where we have used the L^-contraction property of the equation for the last esti- 
mate. Multiplying (|2.ip by a nonnegative test function ip £ {{supp-p, fY) and 
integrating in space and time we have 

u{t)^^ I I {J *{u{t)-l)+)pi^cot f if. 

Jo JR« JR" 

Taking to ■— I/cq, we get Jg^Niuit) ~ 1)9? ^ for aU t e [0,to]. Using an approx- 
imation ipxn where x„ sign^(u — 1), we deduce that J^fi{u{t) — l)+ip — 0, so 
that 

(3.1) suppp,(w(t)) C suppp,(/), t e [0,to]. 



Estimate of the support of u. Lemma 13711 implies then that 
suppp,(wf(t)) C suppp,(/) +Brj, t e [0,to]. 
This means that for any ip G ((supp25/(/) + Bb,,Y) we have 

u{t)f}= f f ut{t)ip = 0, te [0,to], 

Jo JR" 

that is, 

(3.2) suppp,(u(t)) C suppp,(/) + t E [0,to]. 

Iteration. Notice that to depends on the initial data / only through its L^ norm. 
Hence, since the L^ norm of the enthalpy is time invariant, the arguments can be 
iterated to obtain the result for all times. □ 
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Remark. The same argument can be used for initial data which are not compactly 
supported, to show that some positive time will pass before the temperature be- 
comes positive at any given point in the complement of the support of (/ — 1)+. 

The last two results have counterparts for BC-solutions. 

Theorem 3.3. Let f G BC+(M^), and let u he the corresponding BC-solution. 
Then: 

(i) Ut{x, t) = for any x G (supp(i;(-, i)) + BrjY, t ^ 0. 

(ii) //sup|^|^^ /(a;) < 1 for some R> 0, then v(-,t) is compactly supported for 
all t > 0. If moreover f £ Cc(IR.^), then u{-,t) is also compactly supported 
for all t > 0. 

Proof, (i) The proof is similar (though even easier, since the supports are under- 
stood in the classical sense) to the one for L^-solutions. 

(ii) Using Lemma [2. Ill we get 

{J*{u - l) + ){x,t) II J||Li(RiV)||(u - l)+||Loo(RiV) < ||/||l-(R«). 

Therefore, from (|2.ip we have 

(3.3) u(a.,i) ^/(2:)-Ht||/||L»(M«) sup /(x)-|-t||/||Loo(K«), \x\^R. 

Thus, for all |a;| ^ R and t ^ (1 — sup|^|^fl /(a;))/(2||/||Loc(]giV)) we have u{x,t) < 1, 
and hence v{x,t) — 0. Then, by (i), u{x,t) = f{x) for all \x\ ^ R + Rj and 
t = (1 — sup|j.|j,^/(a;))/(2||/||Loo(R]v-)). We now proceed by iteration. □ 

3.2. Equivalence of formulations. As a corollary of the control of the supports, 
we will prove that if the initial data are in Li'^(IR^) n Co(K-^), then the L-'^-solution 
is in fact continuous. 

Proposition 3.4. Let f £ L^(R^) n Co(]R^). The corresponding -solution is 
continuous. 

Proof. We start by considering the case where / is continuous, nonnegative and 
compactly supported. Since a BC-solution with a continuous and compactly sup- 
ported initial data stays compactly supported in space for all times, it is also inte- 
grable in space for all times. Moreover, u € C([0, T]; L^(K^)) . Hence, it coincides 
with the L^-solution with the same initial data. 

We now turn to the general case, that will be dealt with by approximation: let 
{/«} be a sequence of continuous and compactly supported functions such that 

ll/ri - /IIl«=(K") < Win - /||l1(R«) < 

Let M^, be the L^-solutions with initial data respectively /„ and /, and u^, 
u'^ the corresponding BC-solutions. We know that — u^. Now, using the L^- 
contraction property for L^'^-solutions, we have that ||w^ — u^||li(k"x[o,t]) ~^ for 
any T e [0,oo). Moreover, by Lemma B. 101 ||u^ — u'^||L~(RiVx[o,T]) — ^ 0- Hence the 
result. □ 
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We can now use an approximation argument to prove a comparison principle for 
general initial data in BC+(M^). This can in turn be used to show that initial data 
in that class yield one-phase solutions (satisfying u ^ 0). 

Corollary 3.5. Let fi,f2 G BC+(]R^), andui, U2 the corresponding JiC- solutions. 
If h sS h then ui s^U2. 

Proof. Let {/2,n} be sequences of nonnegative, continuous and compactly 

supported functions such that /i.„ /i, /„_2 — >■ /2 uniformly, and ^ f2,n- 
Since /i,„, f2.n G L^(M^) n Co(]R^), the corresponding BC-solutions U2,n, are 
also L^-solutions. Hence, the comparison principle for L^-solutions, Corollary 12. 5[ 
yields ^ U2,n. Passing to the limit in EC, and using the continuous dependence 
of BC-solutions on the initial data. Lemma r2.10[ we conclude that ui ^ U2- □ 

3.3. Retention for u and v. We next prove that the supports of both u and v 
are nondecreasing. We denote this property as retention. 

We start by considering the case of BC-solutions. 

Proposition 3.6. Let f G BC+(M^), and let u be the BC-solution to problem (|l.ip . 
Then dfU ^ — u and dtv ^ —v for all t ^ 0. Ln particular, u and v have the retention 
property. 

Proof. We have 

dtu ^ J * {u — 1)+ — (m — 1)+ ^ —u, 
which, after integration, yields 

(3.4) u{x,t) ^ u{x,s)e~^^-''\ t ^ s. 

This implies retention for u. 

Concerning v — {u — 1) + , we have 

dt{u - 1)+ = dtU ■ t{u>i} > -~{u - 1) + , 
that is, dtv ^ — u, from where retention follows. □ 

For L^-solutions we also have retention for both u and v. In this case the supports 
have to be understood in the distributional sense. 

Proposition 3.7. Let f G L^'^(M^), and let u be the -solution to problem (jl.ip . 
Then u and the corresponding v have the retention property. 

Proof. Let {/„} be a sequence of functions in C^(K.^) such that ||/„ — /||li(r«) ^ 
0. Let Un be the L^-solution (which coincides with the BC-solution) with initial 
data /„. Thanks to the L^-contraction property, ||un(^) — '^'(^)IIli(r~) ~^ for all 
t > 0. Moreover, since the temperature w is a Lipschitz function of u, we also have 

\\Vn{t) - w(t)||Li(R«) 0. 

Let if be any nonnegative function, G C(?°(R^). Multiplying (|3.4p (with u„ 
instead of u) by (^9, integrating and letting n — > oo, we get 

/ u{t)ip ^ e / u{s)ip, t > s, 

JR« JR" 
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from where retention for u in the distributional sense is immediate. The argument 
for V is identical. □ 

3.4. Localization of the supports. Our next aim is to prove, in the case of a 
nonincrcasing kernel J, that for a wide class of initial data the supports of both u 
and V are localized: they are contained in a ball of fixed radius for all times. 

The result will follow from comparison with solutions with initial data in Co(K''^) 
that are radial and strictly decreasing in the radial variable. Such solutions are 
continuous and radial (the latter fact comes from the uniqueness of the solutions 
and invariancy under rotations of the equation). Moreover, for any time t ^ 0, 
the water zone {v{t) > 0} is compactly supported (see Theorem 13. 3p . The main 
technical difficulty stems from the fact that we are not able to prove that these 
solutions are decreasing in the radial variable for all times. 

Lemma 3.8. Let J be nonincreasing in the radial variable and f G Co(]R^) non- 
negative, radial, and strictly decreasing in the radial variable. Then the support of 
v{t) is a ball of radius r(t) for every t ^ and the function r is continuous on 
[0,oo). 

Proof For t > 0, let r{t) := inf{r > : supp(w(-,i)) C 5^(4)}. Thanks to The- 
orem 13.31 this quantity is well defined. By the retention property for v (Proposi- 
tion l3.7p . the function r is nondecreasing. Notice that a priori supp(?;(t)) could be 
strictly contained in though we will prove that this is not the case. 

Continuity of r. Assume for contradiction that r(iQ ) < ^(ij) at some time 
to > 0. For any x such that \x\ > r^f^) we have v{x,t) = for all t ^ to- If 
moreover |a;| = '''{to)' continuity of u yields u{x,to) = 1. 

Let Xa = (a, 0,...,0), x^ = (6, 0, ...,0), with r{to) < a < b = rif^). We 
consider w{t) :— u{xa,t) — u{xf,,t). For any t to, we have 

w'{t) = ((J * v){Xa,t) - v{Xa,t)) - ((J * v){xb, t) - v{xt, t)) 
= ( J * v){Xa,t) - (J * v){xt,t) 

= _ v{y,t){j{xa-y) - J{xb~y))Ay. 

•'\v\<r(to ) 

Since \xa — y\ < \xf, — y\, for all \y\ ^ ''(^o ), then J{xa ~ y) ^ J{xb — y) in this 
region, because J is radially noincreasing. Thus we obtain w' ^ Q. Using that 
u>(0) = f{xa) — f{xb) > 0, we obtain w{to) > 0, so that u(xa,io) > u{xb,tQ) — 1. 
This is a contradiction, since v{xa,to) = 0. 

Continuity at i = is easier. On one hand, from the retention property we 
have r(0+) ^ »'(0). On the other hand, since / is strictly decreasing, we have 
f{x) < 1 if and only if \x\ > r(0). But, by the continuity of u, f{x) = 1 if 
a;| — r(0+). Therefore, r{Q~^) cannot be strictly greater that r(0). We end up with 
r(0) = r(0+). 

CONECTEDNESS. We now prove that supp(i;(-, t)) is connected for all posi- 
tive times. Assume, on the contrary, that there exists some > such that 
supp(u(-, i*)) is disconnected. Hence, there are values r(0) < a < b < r{t^) such 
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that v{x,t^:) = if a ^ |x| < 6 and v{x,t^,) > if 6 < < b + S for some 5 > 0. 
The retention property for v imphes that v{x, t) = for a ^ |a::| ^ 6, ^ t < t*. 

Let td G (0,t*) be the time when the disconnected region outside the bah Bi, 
appears, 

td ■■= sup{t > : v{x,t) = for \x\ > 6} = sup{t > : v{x,t) = for |a:| ^ a}. 

Obviously, r{td) ^ a and, on the other hand, r(tj) ^ b > a, a contradiction with 
the continuity of □ 

Lemma 3.9. Let J be nonincreasing in the radial variable. If ^ f ^ g a.e. for 

some g £ L|'^(M"'^) H Co(M^) radial and strictly decreasing in the radial variable, 
then there exists some R, depending only on g, such that supp(?j(i)) C Br for all 
t ^ 0. 

Proof. Let Ug be the L^-solution with initial datum g and Vg = {ug — 1)+. By com- 
parison, supp(u(i)) C supp(wg(i)). Lemma [3.81 implies that supp(wg(t)) = Brg(t)- 
The radius rg{t) can be estimated using the conservation of mass, 

5 = / Ug{t) ^ / 1 = I SUpp(Wg(t))| = tJjVrg(t)^, 

"'R" J{Ug{t)>l} 

where wjv is the volume of the unit sphere in R^. This implies that 

SUpp(v(t)) C SUpp(Vg(t)) C Br, i? = / g/uJN I 

\JR" / 

□ 

As a corollary of the localization of the support of the temperature, we obtain 
that ||w(i)||Li(jjjv) tends to zero as i — ^ oo with an exponential rate. 



Corollary 3.10. Let J and f satisfy the hypotheses of Lemma l3.9l Then there are 
constants C,k > such that \\v{t)\\i^imN\ ^ Ce~^* for all t^O. 



Proof. Let R be such that supp(w(t)) C Br for all times, and V the L^-solution to 
the nonlocal heat equation in Br, 

Vtix,t) = J*V{x,t) - V{x,t), xe BRt>0, 
V{x,t) =0, x^Br, t>0, 

V{x,0)^if-l)+ix), xeBR. 

Since / is bounded, (/ - 1)+ E 1?{Br). Hence fT2| Theorem 2], ||V"(t)||L2(i3„) 
decays exponentially in time. 

As V is subcaloric. Lemma 12.61 u ^ y in Br x (0, oo). This implies 

v{t) < (w^i?^)i/2||y(i)||L2(s,) < Ce-^K 

□ 

For general integrable data we are only able to obtain a power-like decay rate. 
Corollary 3.11. Let f G Li(R^). Then ||t;(i)||Li(RiV) = 0(t-^/2). 
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Proof. Since v is subcaloric and nonnegative, it is enough to compare it from above 
with the solution, V, to the non-local heat equation (jA.l|) with the same initial 
data. Thus, using the representation formula (jA.2p for solutions to (jA.l[) (see the 
appendix), we get 

/ v{t)= f s;e-*||(/-l)+||Li(R«)+ / c^(t) *(/-!)+, 

JR" J{v{t)>0} J{v{t)>a} 

where w is the regular part of the fundamental solution to (jA.ip . Then we notice 
that the measure of the support of v{t) is uniformly controlled, 

\{v{t) > 0}| = \{uit) ^ 1}| ^ / u{t) ^ ||/||li(r«). 

Thus, using that ||a;(t) ||loo(rn-) ^ Ct~^^^ (see P^), we obtain 

□ 



If the initial data are bounded and compactly supported, we can obtain quanti- 
tative estimates for the supports of u and v. These estimates are sharp, as can be 
checked by considering indicator initial data. 

Lemma 3.12. Let J be nonincreasing in the radial variable and f nonnegative, 
bounded and compactly supported, contained in the ball of radius R f . Then 

SUpp(w(t)) C Bii„, SUpp(M(i)) C S„,ax{fl^,_R„+flj}, Rv = 1 1 / 1 1 ^^{.^'(Riv ) -R/ • 

Proof. To obtain the estimate for the support of v we use Lemma [3.9l with functions 
gn approximating Isj,^. ||/||l°°(r") from above, and then pass to the limit in n. The 
estimate for the support of u then follows from Lemma 13.11 □ 

Remark. A better result should hold with a radius i?„ depending on the mass of 
the initial data above level one, instead of the L°°-norm of /. 

3.5. Creation of mushy regions. Since the L^-solutions to our problem are not 
necessarily continuous, we need a distributional definition of the mushy region. 

Definition 3.13. The mushy region at time < ^ of a nonnegative (L^- or BC-) 
solution u to p.ip is 

(t ) : = Int ( suppi,, {u{t))\ supp^, ( (it (t) - 1 ) + ) ) . 

Remark. When u is continuous, M.{t) = {0 < u{x,t) < 1}. 

Here comes one of the main features of our model: it allows the creation of mushy 
regions. 

Theorem 3.14. Let f G Li'^(R^) be a nontrivial initial data such that suppj5;(/) = 
supp-p,((/ - 1)+). Then, 

A^(t) = {0<dist(a;,suppp,(/))<i?,7}, ie[0,io], = \rjT, ■ 

II-'IIl°°(r")II/IIli(r") 
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Proof. Wc first observe that the assumptions on the initial data imply that there 
are no mushy regions initially, A^(0) = 0. Moreover, 

(3.5) snppj,,{v{t)) = supp-p,(/), t e [0,to]. 

The upper inclusion is just and the lower one follows from the retention 

property and the equality of the supports of / and (/ — 1)+. 

The inclusion A4{t) C {O < dist (a;, suppp, (/)) < i?,/} is an immediate conse- 
quence of equations ()3.2|) and p.5|) . 

Let us turn then to the other inclusion. Let tp he a nonnegative and nontrivial 
test function compactly supported in {0 < dist(x, suppx),(/)) < Rj}. Using p.ip . 
we have 

/ u{t)^= [ [ [ J*v{t)^~ [ [ v{t)^ 

jRi^ Jr« Jo Jr" Jo Jr" 

= 11 J*v{t)ip, t e [0,to]. 
Jo Jr" 

Since supp(J * v{t)) = suppxi,{{v{t)) + Bn, = snpp-p'if) + ^B.j, we conclude that 
/ u{t)ip= J*v{t)ip>Q, te[0,to]. 

JR« JR« 

In other words, 

{0 < dist(a;,suppp,(/))) < Rj} C supp^,, t G [0,io], 

which combined with (|3.5p gives the required inclusion. □ 



3.6. Emergence of disconnected water regions. Another interesting feature 
of our model is that disconnected components of water may appear suddenly at 
a positive distance from the already existing water components. This is another 
example of a phenomenon that occurs for the non-local model but not for the 
local one. The reason is that, contrary to the local model, problem (|l.ip allows 
middle-range interactions, (up to a distance Rj). 

Let us now construct examples exhibiting this phenomenon. We will keep things 
simple, so that the underlying mechanism is better understood. But the result can 
be easily generalized to more complex situations. 

Let us assume that / is a bounded and continuous initial data with three well 
differentiated zones: 
(3.6) 

a "warm" water zone, W, where the enthalpy is above 1; 
< a low-enthalpy ice zone, II, where / is clearly below 1; 

a "high" -enthalpy ice zone, Ih: where / is close to, but below level 1. 

We shall see that if the enthalpy in Ih is close enough to 1 and this zone is not 
too far from W, then Ih melts before the low-enthalpy zone does. Therefore, if II 
"separates" W and Th initially, a new disconnected component will emerge in the 
water zone. It is enough to prove it for Th — {x}, the general case resulting from 
this. 

Definition 3.15. A set iS C separates A and B if: 
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(i) 5^ has at least two open connected components; 

(ii) A and B lie in two different connected components of 5^. 

Theorem 3.16. Let J > m Br,, x G and consider two non-empty sets 
W,Il C such that 

• W is open; 

• dist(a;, W) < Rj; 

• Tl separates W and {x}. 

Let f e BC+(R^) such that 

• f > I inW, f <l in ; 

• O^/^l — T^inlL for some fixed rj G (0, 1); 

Then there exists e G (0,?]) such that if 1 — e < f{x) < 1, then in a finite time 
there appears in the water zone a new connected water component Cx containing x. 

Figures [Hand [2] illustrate the phenomenon, the exact meaning of notations being 
defined within the proof that follows. 




Figure 1 . High-enthalpy zone near melting. 




Figure 2. New disconnected component in water zone. 



Proof. We proceed in three steps as follows. 

Step 1 - i^or < i < 77/|1/|1l~(k«): II n {v{t) > 0} = 0. 
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Let y £ II and t < f?/||/||L°°(R«)- Using p.3p . we have u{y,t) < 1. Thus, such 
an y remains in the ice zone. 

Step 2 - Given < t < ?7/||/||l°°(r"); there exists Eq > such that for all 
e e (0,eo); we have x G {v{t) > 0}. 

Assume on the contrary that v{x,t) = for all e > 0. By our assumptions, 
since dist(x, {/ > 1}) < Rj, there exists < p < Rj and yo E {f > 1} such that 
dist(a;, yo) < p. Taking 6 = {f{yo) — l)/2 > 0, we have f{yo) > 1 + <5 so that the set 

/(-i)((l + <5,+oo))={/>l + <5} 

is open and contains yQ. Similarly, the ball Bp{x) is also open and contains yo so 
that the intersection {/ > l + 5}r\Bp{x) contains at least a ball i?r(2/o) centered at 
yQ with a positive radius r > 0. This has two consequences that we use below: first, 
{/ > 1 + (5} n Bp{x) has a positive Lebesgue measure; second, x — B^iyo) C -B_r,j, 
so that by assumption on J (up to taking a r' < r), J(x — y) is uniformly bounded 
away from zero on Br{yo)- 

Then, using the retention property for u. Proposition 13.61 we have 

u{x,t)>l-e+ (J * s) ds ^ 1 - e + e"* / {J * v){x,0) ds. 
Jo Jo 

We estimate the integral as follows: 

(J*w)(a;,0) ^ / Jix-y)6dy 
J{f>i+s}nBf,(x) 

^ (5 • min Jix~y) ■ \{f > l + 6}nBp(x)\ = C > 0, 
BAyo) 

with C independent of e. Hence, if e > is small enough we get 

u{x,t) > 1 - e + ie~*C > 1, 

which is a contradiction. 

Step 3 - For i and e as above, there is a new connected component containing x 
in the water zone. 

Since x and W are separated by II , there exist two open sets Oi , O2 such that 
Oi U O2 C [IlY, Oi n ©2 = and 

{x}eOi, WcOs. 

Thus, initially all the water zone is contained in O2 while x belongs to Oi. 

Now, by the retention property, we know that at time i the water zone still 
contains W. But at time t the water zone also contains x. Hence we are in the 
following situation: 

{v{t) > 0} D {x} U W, {v{t) > 0} C [XlY. 

At time t, let us denote by Cx the connected water component containing x and by 
Cyy the one containing W, which are both non-empty. Since {v{t) > 0} C {IlY = 
Oi U O2, the union being disjoint, it follows that necessarily Cx C Oi, Cy^ C ©2. 
Hence we deduce that 

n Cvv = 0. 
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In other words, at time t, a new component has appeared in Oi which was not 
present initially, and it is even disconnected from all the water zones in 02- D 

Remark. A similar phenomenon takes place for non-continuous, integrable data. 
This can be proved either by using L^-theory techniques, or by approximation from 
above and from below with continuous data. 



4. The local Stefan problem as a limit in the macroscopic scale 

If the support of the kernel is shrunk to a point through a suitable rescaling, we 
recover the local model. For the case of a bounded domain with Neumann boundary 
data, such convergence was already considered in [3] in the abstract setting of semi- 
group theory. We will give here an alternative, more direct proof, adapted to our 
problem. In addition, we will prove that mushy regions disappear in the limit. 

Given a fixed initial datum / G L^(R^) n L°°(R^), we consider the problem 

(4.1) dtu' ^\{J,*v' -v'), v' = {u'^l)+, u'{-,0) = f, 

where = J{- / e). Since Jg is a unit mass kernel, compactly supported in 
the ball B^^j, the various properties of solutions of (|l.ip that we derived in the 
previous sections are still valid for solutions of (|4.ip . This latter problem admits 
a weak formulation, which will show to be quite convenient when passing to the 
limit: for any test function e C^(]R^ x [0, oo)) we have 

(4.2) / u%mt)= f fm+ f I [dt^W + ^f f {Js*cf>-cj,X. 
This follows from Fubini's theorem, since is symmetric. 

Lemma 4.1. Let f G L-'^(M^) nL°°(M^). The family {u^} is relatively compact in 
Li,,(]R^x(0,^)). 

Proof. We first prove the relative compactness of {u^(t)} in Lj'Q^(M^) for alH > 
by means of Frechet-Kolmogorov's compactness criterium. To this aim we use the 
L^-contraction property to show that: (i) the functions u^{t) are uniformly bounded 
in L-'^(R^); and (ii) for any compact set oj C M.^ 

[ \u^x + h,t)-u'{x,t)\dx ^ I \f{x + h)- f{x)\dx^Oh{l). 

Juj JR« 

where Oh{^) tends to as — > independently of e. Hence, along a subsequence, 
u^it) u{t) in L[^^(K^), for some function u{t) e Li(M^). We infer that it'^(a;,t) 
converges for almost every {x,t). 

Since ||u'^(i)||oo ^ ll/llcx), we may now use the dominated convergence theorem 
to prove that converges to u in L|qj,(M^ x [0, cxd)). □ 

This compactness result gives convergence along subsequences. The possible 
limit functions turn out to be weak solutions to the local Stefan problem 

(4.3) = !^A(w - 1)+, u{;Q) = f, 
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where m2 := /jjn |zp J(2;) dz is the second-order momentum of the kernel J, which 
is finite, since J is compactly supported. Since this problem has a unique weak 
solution [I], convergence is not restricted to subsequences. 

Theorem 4.2. Let f € Li(R^)nL°°(M^). The sequence {u^} of solutions to (j4l|) 
converges as e ^ in Ll^^{M.^ x [0,oo)) to the unique weak solution of (|4.3p . 

Proof. Along a subsequence, {w^} converges strongly in Lj'Q^(K^ x [0, oo)) to some 
function u, see Lemma [4.11 We also have convergence for {v'^} along some sub- 
sequence to V = {u — 1) + . Since the L-'^-norms of solutions to (|4.1[) do not in- 
crease with time, see Corollary 12. 5[ the sequence {we} is uniformly bounded in 
L°°((0,oo);Li(R^)). 

If we perform a Taylor expansion and use the symmetry of J, we get 

J{z) (^4){x) — 4>{x ~ £z)^ dz 

= :^A</. + o(l) ase^O+ 

uniformly in (recall that is compactly supported and smooth). Hence, 

passing to the limit in (|4.2p . we get that u is the unique weak solution to the local 
problem: for any test-function </) G C^(M^ x [0,oo)), 

/ umt)^ f fcj>{Q)+ f f cl>,u+^ f f iA4>)iu-l)+. 

JR« JR" JO Jr" ^ JO Jr" 

□ 




We next study the limit e — > for the mushy region A4'^{t) associated to for 
initial data such that A^'^(O) = 0. We first estimate the size of the mushy region 
for fixed e > 0. 

Theorem 4.3. Let f G Lil^(R^) be a nontrivial initial data such that suppj5/(/) = 
suppx)'((/ - !)+)■ For any t > 0, 

(4.4) M^t) C {x e : < dist{x,suppj,,{v%- ,t)) < sRj}. 

Proof. Let x ^ suppp/(t;'^(-, t)) + i3e_Rj for some fixed t > 0. Since the support of 
w^(t) is nondecreasing, then for any ^ s ^ t, x ^ suppp,(w'^(-, s)) -I- B^jij. Then 
Lemma [5TT] implies that u'^{x,t) — f{x) = 0, which implies (|4.4[) . □ 

Unfortunately, since the convergence of the functions {m^} is rather weak, this 
is not enough to prove the convergence of the mushy regions. However, we will be 
able to prove that their limsup, 

M*{t) = limsup7W"(t) = n U -^"W' 

^'^^ »j>0))<e 

consisting of all points x such that for any 77 > there exists an e £ (0, 77) such that 
X £ M"^ (t) , is a. negligible set. 

Recall that for the local problem, under our assumptions on the initial data 
no mushy regions are created, so that the supports of v{t) and u(t) coincide for 
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all t ^ 0. Moreover, v is continuous in M.^ x (0, oo) [9], hence the distributional 
support of v{t) can be understood as the closure of the set {v{t) > 0}. 

Corollary 4.4. Let f e L\{R^) such that suppp,(/) = suppp,((/ - 1)+). Then, 
for any t > 0, A4*{t) has zero N -dimensional Lebesgue measure. 

Proof. We know that u'^{-,t) converges pointwise to u{-,t) except on a set F which 
has zero A^-dimensional Lebesgue measure. Apart from the set F there are three 
possibilities: 

(i) X e d{v{t) > 0}. This set is known to have zero A^-dimensional Lebesgue 
measure [TD]. 

(ii) X e {v{t) = 0} \ d{v{t) > 0}. Then since there are no mushy regions for the 
limit (local) equation here, necessarily, u{x,t) = 0. This implies that for e small 
enough (say less than eo), we have u^{x,t) < 1, thus v^{x,t) = 0. Hence x cannot 
belong to M*{t) because x belongs to no mushy region for e < eg, see Theorem l4.3l 

(iii) X e {v{t) > 0}. Then for e small enough, u^{x,t) > 1 because it converges 
to u{x,t) > 1. Thus, for e small enough, such an x does not belong to any mushy 
region, hence it is not in the limsup 

Thus we have proved that M*{t) is included in F U d{v{t) > 0} which is a 
negligible set. □ 



Our next aim is to describe the large time behavior of the solutions to our model. 
For the local Stefan problem it is given by a 'mesa'-type problem [T7]. To be more 
precise, u converges to f — f + Aw, where w solves the elliptic obstacle- type 
problem 



In our case the limit is also given by a 'mesa', but now of a non-local character, 
see below. This is to be contrasted with the large time behavior of the non-local 
heat equation in the whole space, which is given by the solution of the local heat 
equation with the same data, and hence by a multiple of the fundamental solution 
of the latter equation. 

5.1. Formulation of the Stefan problem as a parabolic non-local obstacle 
problem (in complementarity form) . We consider here a nonnegative initial u 
giving rise to a nonnegative solution u. To identify the asymptotic limit for u, we 
define the Baiocchi variable 



A variable of this kind was first used by Baiocchi in 1971 to deal with the dam 
problem [4] , [5] . The enthalpy and the temperature can be recovered from w through 
the formulas 



5. Asymptotic behavior 



w^O, 0^/ + Aw^l, {f + Aw-l)w = 0. 




(5.1) 



f + J * w — w. 



V = dtw. 
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where the time derivative has to be understood in the sense of distributions. More- 
over, 

^ u — V ^ 1, {u — 1 — v)v = a.e. 
The distributional supports of v and w coincide for all times. 
Lemma 5.1. For any t > 0, we have swppj^,{v{t)) = suppj)/(z«(t)). 

Proof. Let x G supp-p, {w{t)). Then, given ip G C^(R^), ip ^ 0, (p{x) > 0, we have 

< / w{t)ip ^ f ( f v{s)ip] ds. 

Hence, there exists s < t such that J^^j v{s)(p > 0, i.e., x £ suppj5/(u(s)). Using the 
retention property for v we finally get that x e suppp/ 

Conversely, assume that x e suppj,/ Then, given ip G C^(R^), (p ^ 0, 
(p{x) > 0, we have J^^ v{t)'f > 0. Since v e C([0, oo); L-^(R^)), we have that there 
exists a value ^ > such that Jjjj„ v{s)(p > for all s £ {t — 6,t). This implies that 
/ijN w{t)ip > 0, hence x S suppj,/ (u'(t)). □ 

Thanks to this lemma, (m — 1 — v)w — a.e. Hence w solves a.e. the comple- 
mentarity problem 

(5.2) w ^ 0, 0^f + J*w — w — dfW ^1, {f + J*w — w — \ — dtw)w = 0, 

plus the initial condition w{0) — 0. An analogous formulation for the local Stefan 
problem was given in [Tl], see also [IB] . 

5.2. A non-local elliptic obstacle problem. If \\v{t)\\i^i(js^N^ dt < oo, then 
w{t) converges monotonically and in L^(R^) as t — > oo to 

v{s)ds e L\R^). 



Thus, see (|5.ip . u{-,t) converges point-wisely and in L^(R^) to 

f = f + J *Woo ~ Woo- 
Passing to the limit as i — ?> oo in (|5.2p . we get that Woo is a solution with data / to 
the nonlocal obstacle problem: 

{Given a non-negative data / G L^(R^), find a non-negative 
function w G L^(R^) such that 
0^f + J*w — w^l, {f + J If w — w — l)w ^ a.e. 

This non-local obstacle problem has a unique solution. The proof is based on the 
following Liouville type lemma for J-subharmonic functions. 

Lemma 5.2. Let w G L^(R^) such that w ^ 0, w ^ J * w a.e. Then w = a.e. 

Proof. Assume first that w is continuous, and fix e > 0. Since w is integrable, there 
is a radius R such that 

e 



w < 



x\^R ||>^1|l°=(R«) 
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Hence, for |x| ^ R + Rj 

(5.3) ^ (J * u))(a;) ^ II J||l<»(rn) / w s$ || J||loo(rjv) / w^e. 

JBhj{x) J\x\^R 

So, let us assume that for some x G M^, w{x) > e. Then the maximum of w is 
attained at some point x G Bji^ji, and 

maxw = w(x) > e. 

Using that w ^ J * w, we first deduce that w{x) — 'w{x) in Bij j{x) and then, 
spreading this property to all the space by adding each time the support of J, we 
conclude that w = w{x) > e in all M^. But this is a contradiction with (|5.3I) . So, 
we deduce that ^ it; ^ e for any e > 0, hence w = 0. 

If w is not continuous, we consider Wn = w*p„, where p„ is an approximation of 
the identity. The continuous function Wn satisfies all the hypotheses of the lemma, 
hence w„ = 0. Letting n — oo we obtain w = a.e. □ 

We will also need the following non-local version of Kato's inequality (see [12] 
for the local inequality), 

(5.4) {J * w — w)l[^yoy ^ J * w+ — w+ a.e., 
which is trivially valid for any function in Ljqj.(R^). 
Theorem 5.3. Problem (|OPI) has at most one solution. 

Proof. The key point is that solutions to Problem ()0P|) satisfy 

f^f + J*w~w, f & f3{w) a.e., 

where /3 is the sign graph, see [7 for the local case. 

Let Wi, i = 1,2, be two solutions to (jOP[) with initial data /, and fi be the 
corresponding projections. Since /,; £ l3{wi), we have 

^ (/i - /2)l{tui>tu2} = (J * {wi - W2) - {wi - W2))l{^i>t„2} a.e., 
from where we get, using Kato's inequality ()5.4|) . that 

(wi - W2)+ ^ J * {Wi ~ W2)+. 

Therefore, {wi — 102)+ satisfies the hypotheses of Lemma [5.21 We conclude that 
wi ^ W2. Interchanging the roles of wi and z«2, we get the result. □ 

5.3. The mesa problem. At this point we have a precise characterization of 
the large time behavior of solutions to the non-local Stefan problem (|l.ip when- 
ever lk(OllLi(R") is finite. This is the case, for instance, under the hypotheses 
of Lemma [3.91 see Corollarv 13.101 or for general / e L^(R^) if ^ 3, see Corol- 
lary EH 

Theorem 5.4. Let f G L^(R^), and assume in addition, if N = 1,2, the hypothe- 
ses of Lemma 13.91 . If u is the solution to problem (jl.ip and w is the solution of 
problem (jOPp . then u{t) -^f + J*w-w in L^(M^) as t 00. 
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The map V:f^f = f + J*w~w projects the data / onto a 'mesa'- 
type profile: O^/^Jl, / = lin the non- coincidence set {w > 0}. Notice that 
||/||li(r«) — II/IIli(m")- However, m contrast with the local problem, the projection 
/ is not necessarily equal to / in the coincidence set {w = 0}. 

Up to now we have been able to prove the existence of a solution of (jOPp for 
any / e L?^(M^) only if ^ 3. For low dimensions, TV = 1,2, we have needed 
to add the hypotheses of Lemma 13.91 Hence, for low dimensions the projection 
operator V is in principle only defined under these extra assumptions. However, V 
is continuous, in the L-'^-norm, in the subset of Li||_(]R^) of functions satisfying the 
hypotheses of Theorem 15.41 

Corollary 5.5. Let fi, i = 1,2, satisfying the hypotheses of Theorem \bA\ Then 

ll/l - /2||l1(R«) ^ ll/l - /2||l1(R«)- 

Proof. Since (lOPp has uniqueness, any solution with initial data satisfying the 
hypotheses of Theorem 15.41 can be obtained as the limit as f — oo of the solution 
with the same initial data of the non-local Stefan problem. Hence the result follows 
just passing to the limit as i ^ oo in the contraction property for this latter 
problem. □ 

Since the class of functions satisfying the hypotheses of Theorem 15.41 is dense in 
Li'^(M^), we can extend V by continuity to the whole of this bigger space. Thus, for 
any / G L?|_(M^), Vf is the limit in of {Vfn}, where {/„} is any sequence of non- 
negative, measurable, bounded, and compactly supported functions approximating 
/ in Li(R^). 

Let us notice that, though for any sequence of functions {/«} converging to / in 
L-'^(R^) we have convergence of {Vfn}, we are not able to prove the convergence 
of the corresponding solutions {w„} to a solution of (|0P[) . except under the hy- 
potheses of Theorem 15.41 The main obstacle to prove this convergence is the lack 
of compactness of the inverse of £ (recall C = J * v — v). 

5.4. Asymptotic limit for general data. A simple argument now leads to the 
following characterization of the asymptotic limit of the non-local Stefan problem 
for general integrable initial data. 

Theorem 5.6. Let f € Li'|_(R''^), and u the corresponding solution to problem (jl.ip . 
Let V f he the projection of f onto a non-local mesa. Then t) -> Vf in L^(M^) 
as t — i- CX3. 

Proof. Given /, let {/„} C L^(M^) be a sequence of functions satisfying the hy- 
potheses of Theorem 15.41 which approximate / in L^(R^). Let u„ be the corre- 
sponding solutions to the non-local Stefan problem. We have, 

\\u{t)-Vf\\i^l(RN) ||u(i)-U„(i)||Li(B«) + |iMn(i)-^/n||Li(M«) + ll^/"-^/llLi(K«). 

Using the contraction property for the non-local Stefan problem, and the large time 
behavior for bounded and compactly supported initial data, 

limsup \\u{t) - P/||li(R«) 11/ - /n||Li(R") + 11^/" - ^/IIli(r")- 
t—^oo 

Letting n — > cx^ we get the result. □ 
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6. Numerical experiments 

In order to illustrate sonie of the previous results concerning solutions of , we 
show some numerical experiments. We will take / £ Li'j_(R^) compactly supported, 
and J{x) — 0.75(1 — a;^)+. The space discretization is implemented using the 
trapezoidal rule. For the integration in time we have used an ODE integrator 
provided by Matlab®. 

6.1. Creation of mushy regions. In Figure[3]we illustrate the creation of mushy 
regions. As mentioned in the introduction, this phenomenon is absent in the local 
problem; if there are no mushy regions initially, this is also true for any later 
time. On the contrary, in our non-local model, regions with u between and 1 
appear in the neighborhood of the water region as time passes. This is one of 
the main qualitative features of the model. For this simulation we have taken 
/(x) = 2.1[_i,i]. 
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qI ^ , 1 , , , 1 , 1 

-2.5 -2 -1.5 -1 -0.5 0.5 1 1.5 2 2.5 



Figure 3. Creation of mushy regions. 



6.2. Appearance of disconnected water regions. In the local case, it never 
happens that a new water region appears disconnected from the ones that were 
already present immediately before. On the contrary, this indeed happens some- 
times in our non-local model, as shown in Section 13.61 We have exemplified this 
fact in figures |4] and EJ which correspond to an initial datum which is the sum of 
two characteristic functions, /(x) — 2.5 • l[-i.25.-o.5] + 0.99 • l[o.5,i]- 
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Figure 4. Appearance of disconnected water regions. 



u = 1 

Initial Data 
Time t = 1 




0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 
Figure 5. Appearance of disconnected water regions. Zoom. 



6.3. Behavior as e — > 0. In Figure[6]we illustrate the effect described in Section|4l 
When £ — >■ 0, the solution of the non-local problem converges to the solution of the 
local Stefan problem with the same initial data. Moreover, the mushy regions that 
were created because of the non-local effect disappear. The initial datum is again 
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the characteristic function f{x) — 2 • ij. We compute the solution for e = 1, 
0.5 and 0.2. 



E = 1 

e = 0.5 
e = 0.2 

Local solution 




Figure 6. Convergence to the solution of the local problem as £ — ^ 0. 



6.4. Behavior as t 

of the solutions as t 
initial datum 



oo. Finally we show an example of the asymptotic behavior 
■ oo, which is described in Section [5l We have taken here as 



0, 

(sin(5x))+, 
sin(2a;) + 3, 
0, 

0.3, 



X < -4, 
-4 ^ X 
|x| < 1.5, 
1.5 ^ a; ^ 6, 
6 < X < 6.5, 



1.5, 



so that the convergence to a non-local mesa is more evident. 



7. Conclusion 

We have proposed a non-local model to describe the evolution of a mixture of 
ice and water in an intermediate mesoscopic scale, and derived some properties for 
the solutions which are interesting from the physical point of view, in particular the 
creation of mushy regions. We have also proved that the local Stefan problem is 
obtained in the macroscopic limit, and that mushy regions disappear in this latter 
scale if they were not present initially. 

Theorem 13. 141 does not only assert that mushy regions are indeed created. It has 
another interesting consequence reflected in (|3.5I) : there is a waiting time to = (/) 
until which the support oi v — {u ~ 1)+ remains identical to that of (/ — 1)+. In 
other words, the support of the water phase does not evolve until time tg, though 
the temperature in this phase is decreasing because energy is consumed to break 
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Figure 7. Convergence to the mesa as i — > oo. 

the ice. This waiting time can be interpreted as the time needed to break the 
nearby ice phase in order to convert it into water. We are facing a typical "latent 
heat" phenomenon which is not included in the usual local model. In the case of 
equation (|1.3I) . the water phase begins to move instantaneously and ice is melted 
without any waiting time. 

Summarizing, this nonlocal model seems meaningful and allows several natural 
phenomena at an intermediate scale which are not present in the local model. 

Appendix 

In the spirit of [T2], we prove asymptotic estimates concerning the decay of solutions 
of the nonlocal heat equation 

(A.l) ut = J*u-u, u{0) = f. 

We only assume here that the kernel J is symmetric, nonnegative, with unit mass and a 
finite second order momentum denoted by m2. We prove that the solution u is asymp- 
totically similar to the solution of the local heat equation with the same initial data. The 
gain with respect to the paper [T^] is that our results are valid for general initial data 
/ G L^(R'^), without assuming anything about their Fourier transform or their L°°-norm. 

Notice that the solutions will not eventually enter in the class of data considered in [12] 
(unless they were already there initially). Indeed, they can be written as 

(A.2) u{t) = e-\f + {uj{t) * f) 

with Lj{t) smooth, integrable and bounded, [8]. Hence, u is not bounded ii f ^ L°°(R^). 
However, by subtracting e~*/ to u{t) we ensure that u{t) — e~'/ is bounded, since ijj{t) 
is bounded and / G L^(IR^). Moreover, 

- e"*/llLl(E") ^ ll'^WIlLl(R«)||/||Loo(iiiV) SC IliiWIlLlCHiVjll/llLlCH") < OO, 
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SO that the Fourier transform of the difference is in L^(R^). This is important to go back 
to the original variables after making the computations in Fourier variables. 

Theorem A.l. Let f G L^(R'^) and u he the solution of (|A.1|) with u{Q) = /. Let h he 
the solution of 

h, = ^A/i, ft(0) = u{Q) = /. 
Then, as t ^ oo, there exists a function e{t) — >■ (depending only on J and N) such that 
(A.3) t""^^ max \uit) - e - h{t}\ II/IIl1(m")£W- 

Proof. Following [12j Theorem 2.2], we start from 

n(C,t)=e(^«'-^'*/(C) and MC, i) = e '/(C), withc=^, 
but we make a different estimate using u{t) — e ~*/ instead of u{t), 

l--e-7-/l(e,t)dC = / |(e*(^«'-^)-e-*-e-=l«l^*)/(0ldC 



+ 



e|<r(t) 



(e*(-^(«-i) -e-*-e-^l«l'*)/(C)| d^ 
(e*(-^«)-i)-e-*-e-^l«l'*)/(0| d^ 



= 7 + //. 
As in [12], we split 7 into two parts 

I^f |e-=l«l'*^o(0| dC+ / |e*(-^«)-^' -e-*| |uo(C)|dC = 7i+72. 

The only important modification with respect to the proof of [121 Theorem 2.2] concerns 
the term 72. In [12] this term is estimated by using the L^-norm of /, which is something 
we do not want to do here. 

In order to estimate 72 we use that that J verifies 

J{^) < 1 - c|CI^ + ICI^^(C), with h bounded, h{0 as ^ -> 0. 

Hence there exist a, D, 5 > such that 

J(0 ^ 1 - D\^f, for Id «; a and | J(C)| 1 - 5, for \^\ ^ a. 

We decompose 72 by considering separately the sets {r{t) 5C 1^ ^ a} and {1^ ^ a}. The 
integration over {\^\ a} is estimated here taking into account the term e~*/, 

/ |e*'^(«-^' -e-j |/(C)|dC = e-* / |e'-^«' -ij |/(C)|dC. 
Using that in this set | J| ^ 1 — 5, we get that 

|e"<«'-l| 1/(01 ^ 1/(01 E^^^^ 

< ^1^(011/(01 E^!;^ 



11=0 



< ie(^-^'V(01ll/llc 
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Hence 

where £2(i) — exponentially fast, and independently of the initial data. 

Summing up (see [121 Theorem 2.2] for the estimates of the other terms of / + //) 

t'^'''||«(t)-e-*/-ftW||Li(K«) s; ||.f||Li(M«)eW^O 

as f oo for some function e which only depends on J and N. This implies (IA.3[l by 
going back to the original variables. □ 
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